WHAT IS CLAIMED: 



1 . A method for correcting errors in at least one of a plurality of pairs of profiles 
{Am, Cm}, Am being an experiment profile, Cm being a reference profile, where m = 1, 2, 
. . M, M is the number of pairs of profiles, said method comprising 

(a) calculating an average reference profile C of reference profiles {Cm}, m = 1, 2, 

...,M; 

(b) determining for at least one profile pair m e {1, 2, . . M} a differential reference 
profile of Cm and C ; and 

(c) generating for said at least one profile pair m an error-adjusted experiment 
profile A m by a method comprising adjusting said experimental profile Am using said 
differential reference profile determined for said profile pair m, thereby correcting errors in 
said at least one of said plurality of pairs of profiles; 

wherein for each m e {1, 2, ., M}, said error-adjusted experiment profile A m comprises 
data set {A'm(k)}, said experiment profile Am comprises data set {Am(k)}, said reference 
profile Cm comprises data set {Cm(k)}, and said average reference profile C comprises data 
set { C (k)}, wherein said data set {Am (k)} comprises measurements of a plurality of 
different cellular constituents measured in a sample having been subject to condition Am, 
said data set {Cm(k)} comprises measurements of said plurality of different cellular 
constituents measured in a sample having been subject to condition C, and wherein k = 1, 2, 

N is an index of measurements of cellular constituents, N being the total number of 
measurements. 

2. The method of claim 1, wherein said steps (b) and (c) are performed for each 
profile pair m. 

3. The method of claim 2, wherein experiment profile Am and reference profile Cm 
are measured in the same experimental reaction. 

4. The method of claim 3, wherein said C (k) is calculated according to equation 

1 M 

wherein said differential reference profile is calculated according to equation 
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C,,Am,k) = C„{,k)-C{k) 



and wherein said error-adjusted profile is calculated according to equation 



A„{k) = A„-C,,Am,k). 



5. The method of claim 4, further comprising 

(d) calculating for each profile pair m an error-corrected experiment profile A"m 
comprising data set {A'm(k)} by combining said error-adjusted experiment profile A'm with 
said experiment profile Am using a weighing factor {w(k)}, k = 1, 2, N, wherein w{k) is 
a weighing factor for the k' th measurement. 

6. The method of claim 5, wherein said error-corrected experimental profile A"m is 
calculated according to equation 



7. The method of claim 6, wherein said weighing factor w(k) is determining 
according to equation 



where avg_bkgstd is an average background standard error. 

8. The method of claim 7, wherein said avgjbkgstd is determined according to 
equation 

\ ^ ( \ ^ ^ 
avg_bkgstd = — Y —Ybkgstd{m,k) 

Nt:i^Mt;:l J 

where bkgstd (m, k) is background standard error of Cm(k). 

9. The method of claim 4, further comprising determining errors { cr„ } of said error- 
adjusted experiment profiles {Am}. 

10. The method of claim 9, wherein said errors are determined according to equation 



{k) = (1 - w{k)) . (k) + w(k) . (k) . 




o-„ (k) = (k) + mixed _al(k)~2' Cor(k) • C7„ (k) • mixed _ (k) 
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where (A:) is the standard error of Am(k), mixed _ {k) is determined according to 
equation 

. o-,(A:) + (M-l)a,,^(/r) 

mixed _ {k) = 

M 



1 • _ 

where C7,,^(A:) = /— £(C„(^)-CW) 



and where Cor(k) is a correlation coefficient between experiment profile and reference 
profile. 

11. The method of claim 10, wherein said Cor(k) is determined according to 
equation 



( ( c(*) v^ 



Cor(k) = CorMax 



-o.s( 
1-e 



av% _ bkgstd 



where CorMax is a number between 0 and 1 . 

12. The method of claim 7, further comprising determining errors { <t^ } of said 
error-corrected experiment profile {A"m}. 

13. The method of claim 12, wherein said errors are determined according to 
equation 

(k) = V[l - Hk)] • (k) + w(k) ^ o-„ (k) 
where (k) is the standard error of Am(k), cr^ (k) is determined according to equation 

{k) = V (k) + mixed _ (A:) - 2 • Cor(k) - (k) - mixed (k) 

wherie mixed _ cr^ (k) is determined according to equation 

. cr^(A:)4^(M^l).a,,^(A:) 
mixed <j„ (k) = 
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1 ^ — 
where a,,^ (k) = /— £ (C„ (k) - C (k)) 



and where Cor(k) is a correlation coefficient. 

14. The method of claim 13, wherein said Cor(k) is determined according to 
equation 



Cor{k) = CorMax 



avg bkgstd 



where CorMax is a number between 0 and 1 . 

15. The method of claim 2, wherein each said pair of profiles Am and Cm is 
measured in a two-channel microarray experiment. 

16. The method of claim 15, wherein said reference profiles {Cm}, m = 1, 2, M, 
10 are measured with samples labeled with a same label. 

1 7. The method of claim 2, wherein at least one of said plurahty of pairs of profiles 
{Am, Cm} is a virtual profile. 

18. The method of claim 14, wherein said plurality of pairs of profiles {Am, Cm} are 
transformed profiles comprising transformed measurements. 

15 19. The method of claim 18, wherein said transformed measurements are obtained 

according to equations 

A„ (k) = fix) = -^^ 1 + d, 

a 

forXA„,(k)>0 

and 



In 

20 C„ ik) = fix) = 1 + d, 



a 

forXC„(k)>0 
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where experiment profile XAn, comprises measured data set {XAm(k)}, said reference 
profile XCm comprises measured data set {XCm(k)}, where d is described by equation 



-In 



d = 



a 



and where a is the fractional error coefficient of said experiment, h is the Poisson error 
coefficient of said experiment, and c is the standard deviation of background noise of said 
experiment. 

20. The method of claim 14, wherein said experiment profile Am and reference 
profile Cm comprises measurements fi:'om which nonlinearity is removed. 

21. The method of claim 20, wherein said measurements from which nonlinearity is 
removed are obtained by a method comprising 

(i) determining an average profile of all experiment profiles {Am} and reference 
profiles {Cm}; and 

(ii) adjusting each Am or Cm based on a difference between said Am or Cm and said 
average profile. 

22. The method of claim 21, wherein said difference is determined using a subset of 
measurements in the profiles. 

23. The method of claim 22, wherein said subset of measurements in the profiles 
consists of measurements that are ranked similarly between an experiment or reference 
profile and said average profile. 

24. The method of claim 23, wherein said comparing in said step (ii) is carried out 
by a method comprising: 

(111) birming measurements in said subset into a plurality of bins, each said bin 
consisting of measurements having a value in a given range; 

(112) calculating mean difference between said Am or Cm and the average profile in 
each bin; 
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(113) determining a curve of said mean difference as a function of values of 
measurements for said Am or Cm> nonlinear Am or nonlinearjCm-^ respectively; 

(114) adjusting Am or Cm according to equations 
^rik) = A^{k)- nonlinear _A^{k) 

5 or 

Cr (k) = C„ (k) - nonlinear _ (k) 
where k = 1, N. 

25. The method of claim 14, wherein each said experiment profile Am and reference 
profile Cm is a normalized profile. 

10 26. The method of claim 25, wherein said normalized profile is obtained by a 

method comprising normalizing experiment profile Am and reference profile Cm according 
to equation 

NA^(k)= '^-^^'^^ 

A.^ 

and 

15 A^C.W= ^'" — '^^ 

tn 

where A^ is an average of profile {Am(k)}, and C„ is an average of profile {Cm(k)}; 
wherein AC is an average of all profiles calculated according to equation 

J M 

27. The method of claim 26, fiirther comprising normalizing errors of said 
20 experiment profile Am and reference profile Cm according to equation 



and 



o-r(A:)=-^ 
A 



— c" — 
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where cr^ (k) and erf (k) are the standard error of Am(k) and Cm(k), respectively, and 
a^'^ik) and al^^(k) are normaHzed standard error of NAm(k) and NCm(k), respectively. 

28. The method of claim 27, further comprising normalizing background errors of 
said experiment profile Am and reference profile Cm according to equation 



and 



Cm 

where bkgstd^(k) and bkgstd^(k) are the standard background error of Am(k) and Cm(k), 
respectively, and bkgstd^^^ (k) and bkgstd^^ (k) are normalized standard background error 
1 0 of NAm(k) and NCm(k), respectively. 

29. The method of claim 28, wherein said and are an average of 

measurements in profile {Am(k)} and {Cm(k)}, respectively, excluding measurements 
having a value among the highest 10%. 

30. A method of correcting errors in a plurality of pairs of profiles {XAm, XCm}, 

15 XAm being an experiment profile, XCm being a reference profile, where m = 1, 2, . . ., M, M 
is the number of pairs of profiles, said method comprising 

(a) processing said profiles to obtain a plurality of pairs of processed profiles {Am, 
Cm}, Am being a processed experiment profile, Cm being a processed reference profile; 

(b) calculating an average reference profile C of reference profiles {Cm}, rn = 1,2, 

20 ...,M; 

(c) determining for each profile pair m a differential reference profile of Cm and C ; 

and 

(d) generating for each profile pair m an error-adjusted experiment profile A m by a 
method comprising adjusting said experimental profile Am using said differential reference 

25 profile determined for said profile pair m, thereby correcting errors in said plurality of pairs 
of profiles; 
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wherein for each m e {1, 2, . . M}, said error-adjusted experiment profile A'm comprises 
data set (A m(k)}, said processed experiment profile Am comprises data set {Ani(k)}, said 
processed reference profile Cm comprises data set {Cm(k)}, and said average reference 
profile C comprises data set { C (k)}, said experiment profile XAm comprises data set 
5 {XAm(k)}, said reference profile XCm comprises data set {XCm(k)}, wherein said data set 
{XAm (k)} comprises measurements of a plurality of different cellular constituents 
measured in a sample having been subject to condition Am, said data set {XCm(k)} 
comprises measurements of said plurality of different cellular constituents measured in a 
sample having been subject to condition C, and where k = 1, 2, N is an index of 
10 measurements of cellular constituents, N being the total number of measurements. 

31. The method of claim 30, wherein said step (a) comprises normalizing each said 
experiment profile XAm and reference profile XCm- 

32. The method of claim 31, wherein said normalizing is carried out according to 
equation 



15 A„(k) = NA„(k)= ^'"^^^^ 

and 



where NAm and NCm denotes normalized experiment and normalized reference profiles, 



respectively, where is an average of profile {XAm}, and XC^ is an average of profile 



20 {XCm}; wherein XAC is an average of all profiles calculated according to equation 

1 M 

XAC = --Y.^XA„^XC„). 

33. The method of claim 32, further comprising normalizing errors of said 
experiment profile XAm and reference profile XCm according to equation 



25 and 



XA. 



XC„ 
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where (k) and cr^^(A:)are the standard error of XAm(k) and XCm(k), respectively, and 
(y^{k) and erf (A:) are normalized standard error of Am(k) and Cm(k), respectively. 

34. The method of claim 33, further comprising normalizing background errors of 
said experiment profile XAm and reference profile XCm according to equation 

and 



^^^^ m 

where bkgstd^ (k) and bkgstd^^ (k) are the standard background error of XAm(k) and 
XCm(k), respectively, and bkgstd^(k) and bkgstd^ (k)2ire normalized standard background 
10 error of Am(k) and Cm(k), respectively. 



35. The method of claim 33, wherein said XA^ and XC^ are an average of 

measurements in profile {XAm} and {XCm}, respectively, excluding measurements having 
a value among the highest 10%. 

36. The method of claim 35, wherein said step (a) further comprises transforming 
1 5 said normalized profiles to obtain transformed profiles. 

37. The method of claim 36, wherein said transforming is carried out according to 
equations 



1? 



J ^' +^-^"^-^Ar.(k) ^2.^c' .NA„(k) + a' . [NA„ (k)f 

TA„ (k) = fix) = ^ + d, 

a 

torNA„,(k)>0 

20 and 



In 



a 



TC„ (k) = fix) = + d, 

a 

forNC„(k)>0 



77 

NYJD: 1505610.1 



where experiment profile XAm comprises measured data set {XAm(k)}, said reference 
profile XCm comprises measured data set {XCm(k)}, where d is described by equation 



and where a is the fractional error coefficient of said experiment, b is the Poisson error 
coefficient of said experiment, and c is the standard deviation of background noise of said 
experiment. 

38. The method of claim 37, wherein said step (a) further comprises removing 
nonlinearity from each said transformed experiment profile TAm and transformed reference 
profile TCm- 

39. The method of claim 38, wherein said removing nonlinearity is carried out by a 
method comprising 

(al) determining an average transformed profile of all transformed experiment 
profiles {TAm} and transformed reference profiles {TCm} ; and 

(a2) adjusting each TAm or TCm using a difference between said TAm or TCm and 
said average transformed profile. 

40. The method of claim 39, wherein said difference is determined using a subset of 
measurements in said transformed profiles. 

41. The method of claim 40, wherein said subset of measurements in said 
transformed profiles consists of measurements that are ranked similarly between an 
experiment or reference profile and said average profile. 

42. The method of claim 41, wherein said comparing in said step (a2) is carried out 
by a method comprising: 

(a2i) binning measurements in said subset into a plurality of bins, each said bin 
consisting of measurements having a value in a given range; 

(a2ii) calculating mean difference between said Am or Cm and the average profile in 




a 



each bin; 
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(a2iii) determining a curve of said mean difference as a function of values of 
measurements for said TAm or TCm, nonlinear TAm or nonlinear JTCm, respectively; and 

(a2iv) adjusting TAm or TCm according to equations 
TA'J' (k) = TA^ (k) - nonlinear _ TA^ (k) 

5 or 

rcr' (k) = TC„ (k) - nonlinear _ TC„ (k) 
where k = 1, N. 

43. The method of claim 42, wherein experiment profile Am and reference profile Cn 
are measured in the same experimental reaction. 

10 44. The method of claim 43, wherein said C (k) is calculated according to equation 

1 A/ 

wherein said differential reference profile is calculated according to equation 

C,^{m,k) = C„{k)-C{k) 
and wherein said error-adjusted profile is calculated according to equation 
15 A„{k) = A„-C,^{m,k). 

45. The method of claim 44, further comprising 

(d) calculating for each profile pair m an error-corrected experiment profile A"m 
comprising data set {A' m(k)} by combining said error-adjusted experiment profile Am with 
said experiment profile Am using a weighing factor {w(k)}, k = 1, 2, N, wherein w{k) is 
20 a weighing factor for the k' th measurement. 

46. The method of claim 45, wherein said error-corrected experimental profile A' m 
is calculated according to equation 

= (1 - Mik)) ■ A„ (k) + w(k) • A„ (k) . 
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47. The method of claim 46, wherein said weighing factor is determining according 
to equation 



-o.s(. 



C(k) '\ 
avg_bkgstd J 



w{k) 

where avgjbkgstd is an average background standard error. 

5 48. The method of claim 47, wherein said avg bkgstd is determined according to 

equation 



avg _ bkgstd = — X 77 2 t>kgstd (m, k) 



f 1 M 



where bkgstd (m, k) is background standard error of Cm(k). 

49. The method of claim 44, further comprising determining errors { cj^ } of said 
10 error-adjusted experiment profile {A'm}. 

50. The method of claim 49, wherein said errors are determined according to 
equation 

W = ^J al (k) + mixed _ (A:) - 2 • Cor(k) - a, (k) • mixed _ (k) 

where a^(k) is the standard error of Am(k), mixed _ a ^(k) is determined according to 
1 5 equation 

. C7„(k) + (M-\).C7^^(k) 

mixed (k) = 



where <7„,(^) = /-^X(C„(A:)-C(A:))= 



and where Cor(k) is a correlation coefficient between experiment profile Am and reference 
profile Cm- 

20 51. The method of claim 50, wherein said Cor(k) is determined according to 

equation 
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-0.5 I _ ^ 

Cor(k) = CorMax 



-o.5( g'^> r 



where CorMax is a number between 0 and 1 . 

52. The method of claim 51, further comprising determining errors { cr^ } of said 
error-corrected experiment profile {A'm}- 

53. The method of claim 52, wherein said errors are determined according to 
equation 

{k) = ^l[l-wik)]-C7lik) + wik)cT,„ik) 

where ct„ (k) is the standard error of Ani(k), cr„ (k) is determined according to equation 

CT„ (k) = yl<jlik) + mixed _ o"^ (A:) - 2 • Corik) ■ a„ (k) • mixed _ a„ (k) 
where mixed _ cr„ (A;) is determined according to equation 

mixed _ C7„ (A:) = — 

M 



where = l-_|;(C„(A:) -C^)^ 

and where Cor(k) is a correlation coefficient. 

54. The method of claim 53, wherein said Cor(k) is determined according to 
equation 



Cor{k) = CorMax 



-o..{ 



avg _ bkgstd 



where CorMax is a number between 0 and 1. 

55. The method of claim 54, wherein each said pair of profiles XAm and XCm is 
measured in a two-channel microarray experiment. 
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56. The method of claim 55, wherein said reference profiles {XCm}, m = 1, 2, . . 
M, are measured with samples labeled with a same label. 

57. The method of claim 56, wherein at least one of said pair of profiles {XAm, 
XCm} is a virtual profile. 

5 58. A computer system comprising 

a processor, and 

a memory coupled to said processor and encoding one or more programs, 
wherein said one or more programs cause the processor to carry out the method of any one 
of claims 1-57. 

10 59. A computer program product for use in conjunction with a computer having a 

processor and a memory connected to the processor, said computer program product 
comprising a computer readable storage medium having a computer program mechanism 
encoded thereon, wherein said computer program mechanism may be loaded into the 
memory of said computer and cause said computer to carry out the method of any one of 

15 claims 1-57. 

60. A method for generating a differential profile A vs. B fi-om differential profiles 
A vs. Ca and B vs. Cb, comprising calculating said differential profile A vs. B according to 
equation 

20 lratioAB{k) = polarityAC • IratioAC(k) - polarityBC • IratioBC(k) 

where A: = 1, 2, . . ., N, is the index of measurements in a profile, N being the total number of 
measurements; wherein lratioAC(k) = l.o%{A(k) / CA(k)}, if PolarityAC = 1, and lratioAC(k) 
= Log{ CA(k) /A(k)}, if PolarityAC = -1, where A(k), and CA(k) are the k'th measurement 

25 fi-om sample A and Ca, respectively; wherein lratioBC(k) = l.o%{B(k) / CaCk)}, if 

PolarityBC = 1, and IratioACfk) = Log{ Csfk) /B(k)}, ii PolarityBC = A, where B(k), and 
CB(k) are the Ar'th measurement fi-om sample B and Cb, respectively; wherein {A(k)} 
representing measurements of a plurality of different cellular constituents measured in a 
sample having been subject to condition A, {B(k)} representing measurements of said 

30 plurality of different cellular constituents measured in a sample having been subject to 

condition B, and {CA(k)} and {CB(k)} each representing measurements of said plurality of 
different cellular constituents measured in a sample having been subject to condition C. 
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61. The method of claim 60, further comprising calculating an error of differential 
profile A vs. B according to equation 



^IratloAB (^) = V ^LiioAC (^) + ^latioBC (^) ^ 2 • CorMoX ^ CJ i^^,,^^^ (^) " IratioBC (^) 

5 

wherein <JiratioAc(k) and C7iratioBc(k) are errors of lratioAC(k) and IratioBC(k), respectively, 
and wherein CorMax is an estimated maximum correlation coefficient between errors of 
A/C and B/C. 

62. The method of claim 60, wherein A vs. Ca and B vs. Cb are experimentally 
1 0 measured profiles. 

63. The method of claim 60, wherein at least one of A vs. Ca and B vs. Cb is a 
virtual profile. 

64. A computer system comprising 
a processor, and 

15 a memory coupled to said processor and encoding one or more programs, 

wherein said one or more programs cause the processor to carry out the method of any one 
of claims 60-63. 

65. A computer program product for use in conjunction with a computer having a 
processor and a memory connected to the processor, said computer program product 

20 comprising a computer readable storage medium having a computer program mechanism 
encoded thereon, wherein said computer program mechanism may be loaded into the 
memory of said computer and cause said computer to carry out the method of any one of 
claims 60-63. 

66. A method for correcting errors in at least one of a plurality of pairs of profiles 
25 {Am, Cm}, Am bciug an experiment profile. Cm being a reference profile, where m = 1, 2, 

. . M, M is the number of pairs of profiles, said method comprising generating for at least 

one profile pair m g {1, 2, . . M} an error-adjusted experiment profile Am by a method 

comprising adjusting said experimental profile Am using a differential reference profile 

generated using Cm and an average reference profile C determined for said profile pair m, 

30 wherein said average reference profile C is an average of reference profiles {Cm}, m = 1, 2, 

. . ., M; wherein for each m € {1, 2, . . ., M}, said error-adjusted experiment profile Am 

comprises data set {A m(k)}, said experiment profile Am comprises data set {Am(k)}, said 
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reference profile Cm comprises data set {Cm(k)}, and said average reference profile C 
comprises data set { C (k)}, wherein said data set {Am (k)} comprises measurements of a 
plurality of different cellular constituents measured in a sample having been subject to 
condition Am, said data set {Cm(k)} comprises measurements of said plurality of different 
5 cellular constituents measured in a sample having been subject to condition C, and wherein 
k = 1, 2, N is an index of measurements of cellular constituents, N being the total 
number of measurements. 
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